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^ ■ ABSTRACT 

Context. We re-analysed the carbon and oxygen isotopic ratios in tlie atmospheres of the 
two bright K giants Arcturus (a Boo) and Aldebaran {a Tau). 
^ I Aims. These stars are in the evolutionary stage following the first dredge-up (FDU). 

Previous determinations (dating more than 20 years ago) of their '*0/"*0 ratios showed 
a rough agreement with FDU expectations; however, the estimated ^^0/"0 and '"C/'^C 
ratios were lower than in the canonical predictions for red giants. Today these anomalies 
. are interpreted as signs of the occurrence of non-convective mixing episodes. We therefore 

re-investigated this issue in order to verify whether the observed data can be reproduced in 
this hypothesis and if the rather well determined properties of the two stars can help us in 
fixing the uncertain parameters characterizing non-convective mixing and in constraining 
its physical nature. 

Methods. We used high-resolution infrared spectra from the literature to derive the '-C/'^C 
, and '*0/ "O/ "*0 ratios from CO molecular lines near 5 pm, using the LTE spectral synthe- 

■ sis method. We made use of the recently published ACE-FTS atlas of the infrared solar 
5^ . spectrum for constructing an updated atomic and molecular line lists in this spectral range. 

We also reconsidered the determination of the stellar parameters to build the proper atmo- 
spheric and evolutionary models. 

Results. We found that both the C and the O isotopic ratios for the two stars considered 
actually disagree with pure FDU predictions. This reinforces the idea that non-convective 
transport episodes occurred in them. By reproducing the observed elemental and isotopic 
abundances with the help of parametric models for the coupled occurrence of nucleosynthe- 
sis and mass circulation, we derived constraints on the properties of non convective mixing, 
providing information on the so far elusive physics of such phenomena. We find that very 
slow mixing, like that associated to diffusive processes, is incapable of explaining the ob- 
served data, which require a rather fast transport. Circulation mechanisms with speeds 
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Table 1. Oxygen isotopic ratios after FDU at solar metallicity 
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575 


1.25 


1784 


587 


1.3 


1480 


597 


1.4 


1095 


613 



intermediate between tiiose typical of diffusive and of convective mixing sliould be at play. 
We however conclude with a word of caution on the conclusions possible at this stage, as the 
parameters for the mass transport are rather sensitive to the stellar mass and initial com- 
position. At least for a Boo, reducing the uncertainty still remaining on such data would be 
highly desirable. 

Key words, stars: abundances - stars: individual: Arcturus, Aldebaran - stars: late type 



1. Introduction 



Red giant branch (RGB) stars undergo evolutionary stages that start at the so-called First Dredge- 
Up (FDU), a convective mixing process carrying to the surface nuclei from internal layers, previ- 
ously affected by CN cycling. The FDU occurs as the He-core contraction after the Main Sequence 
(MS) is accompanied by a downward envelope extension. It is now well established that this leads 
to a decrease in the '^C/'^C ratio with respect to the MS value (~ 89 in the solar ca se), down to 



Weiss et al 



value s in the range 15-30 (depending on the initial mass and metallicity of the star, see 
20001) . In addition, the carbon abundance drops in the envelope, while that of nitrogen increases. If 
the stellar mass does not exceed ~ 2 M©, the '^O abundance remains unaltered, while that of '**0 is 
mildly reduced. The isotopic ratios "'0/'^0/"*0 expected by the models then lie on a characteristic 
line; their values depend on the initial stellar mass and have m oved strongly in recent years as a 



Palmerini et al 



201 lai especially 



consequence of changes in basic reaction rates (see for example 
their Figure 3). The present situation for these ratios as a function of the stellar m ass, updated with 



2009 



2011 ) and with the last 



2011) is summarized 



the last version of the FRANEC evolutionary code (iCristallo et al 
recommendations available for the relevant reaction rates (lAdelberger et al. 
in Table □ 

This standard picture is challenged by a large amount of abundance determinations (see for 



Brown & Wallerstein 


1989 


Gratton et al. 


2000 



Charbonne] 



2004 



Grundahl et al. 



2002) 



example 

in field and globular cluster low-mass giant stars, showing very low '^C/'^C ratios, sometimes 
almost reaching down to the equilibrium value of the CN cycle (~ 3.5). Anomalies in the C and 
O isotopes were found also in pre-solar C-rich and O-rich g rains of stellar origin, preserved in 



Amari et al. 



2001 



meteorites (for example 
remarkable '^O destruction (iNittler et al. 



Nittler et al.ll2008h . In particular, AI2O3 grains reveal 
I997I) . Some families of this cosmic dust display also 



isotopic shifts in heavier elements including Mg an d Al and others have anoma lies reaching up to 
neutron-capture elements beyond iron (for example]] 



Nicolussi et al.lll997 



1998). 



Observationally, evidence of anomalies is often found in low mass red giants (< 2.3 Mq) 
for phases subsequent to the so-called Bump of the Luminosity Function (BLF), when the ad- 



2 



C. Abia et al.: Carbon and oxygen isotopic ratios in Arcturus and Aldebaran: 



vancing H-burning shell erases th e chemical discontinuity left behind by the first dredge up 



dCharbonnel & Balachandran 



20001). This homogenization facilitates the occurrence of transport 



phenomena; hence, it became common to attribute t he chemical anomalies t o the occurrence of 



episodes of matter cir culation in "conveyor belt s" ( Wasserburg et al 



Palmerini et al 



201 lah or in diffusive processes (IDenissenkov et al 



1995 



1998 



NoUett et al 



Eggleton et al 



2003 



200^ 



These phenomena (going under the names of deep mixing, extra-mixing, or cool bottom processes) 
would link the envelope to regions where proton captures take place, thus accounting for the obser- 
vation that the photospheric material has undergone extensive processing. 

Among the propos ed physical causes for mixing mechanisms one ca n mention rotation itself 



throu gh shear eff ects dZahn 



1992; 



Weiss et al 



tion (ITalon 



20051): gravity wav es (IDenissenkov & Tout 



2007: IDenissenkov et al 



2000; 



Charbonnel 



20041) and meridional circula- 



20031); magnetic buoyancy (IBusso et al 



20091) ; and m olecular weight inver sions leading to heavier materials 



fal ling down in a ligh t er env ironment (lEggleton et al.l l2006h . This last process was identified 
by 



Charbonnel & Zahn 



( l2007b a s the known th ermohaline double-diffusion also occurring in the 



oceans and previously studied bv lUlrichl(ll97ll) in astrophysical environments. 



Although all these physical phenomena may have a role in the complex dynamics hnking red 
giant envelopes to their underlying radiati ve layers, it is not clea r today in which evolutionary 
stage each of them works more efficiently (lUttenthaler et alJl2007h . which is the range of stellar 
masses affected and, therefore, which of them is more suited to explain the abundance changes in 
RGB stars. Even the requirem ent that the mixing episodes occur after the BLF has been recently 



questioned dDrake et al 



201 ih . 



In general, the mentioned elusive processes are not treated in canonical stellar models; in some 
cases they are intrinsically linked to the stellar rotation or to the development of dynamical insta- 
bihties, thus requiring at least two dimensional hydro-codes to be properly modeled. However, the 
use of 2D schemes for general stellar evolution is (at best) in its infancy. 

One has also to notice that clarifying the physics that is behind the chemical pecuharities is 
made difficult by the fact that the interpretation of observations is usually hampered by uncertain- 
ties in fundamental parameters of the chosen stars (stellar mass, luminosity, etc.). The sources 
we consider in this paper, Arcturus (a Boo) and Aldebaran {a Tau), are K-type RGB stars of 
nearly solar mass; they are very bright, are situated at a limited distance from the Sun and are 
well observed. Therefore, they may be less affected than others from this last difficulty, and can be 
considered as good references for studying stellar evolution and spectroscopic abundances in first- 
ascent red giants. Indeed, the determination of their observed parameters (luminosity, radius and 
effective temperature) is rather reliable and their large brightness facilitates the task of getting high- 
resolution, high signal-to-noise-ratio spectra. This is so for both optical and infrared wavelengths, 
thus allowing for an accurate abundance analysis. 

It has been known for decades that the '^C/'^C ratios of a Boo and a Tau share the problems 
discussed a bove for common re d gi ants, being lower than pred icted by the first dredge-up (here- 
after FDU). 



Hinkle et al. 



(1 19761) and 



Tomkin & Lamberj (119841) early derived ratios of 7 and 12 for 



these stars, respectively. These first estimates were then confirmed by subsequent works (see for 
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Peterson et al Thus, the anomaly in their C isotopic ratios 



Smith &Lamber J 19901: 

T 



is a quite robust result 

On the other hand, a previous determination of the oxygen isotopic ratios in these stars 



(IHarris & Lambert 



1984 . using CO lines at 2 and 5 ^im, indicated '^O/'^O- 1100, ""O/'^O- 



550 for Arcturus, and '''O/'^O- 660, "'O/'^O- 475 for Aldebaran. The uncertainty was about 
40 - 50%!. Actually, the mentioned authors noticed the difficulty of explaining simultaneously the 
^^C/'^C and '^O/'^O ratios in these stars within the framework of the canonical FDU models. They 
proposed several possible solutions: among others, strong mass loss prior to the FDU, slow mixing 
during the main sequence and/or a reduction in the rate for ^^0(p, a)'^N. No satisfactory solution 
was however found, so that the problem remained ope n since then. 



Recent theoretical calculations at solar metallicity jPalmerini et al.ll201 lal) . including revisions 
of critical nuclear rates 0, for masses close to 1.2 Mq, found O isotopic ratios shown in Table 10. 
A quick inspection to Table 1 reveals that, considering the observational uncertainty in the Harris 
& Lambert data for oxygen, the '^O/'^O ratios in both stars can be considered to be in rough 
agreement with the new theoretical predictions of stars with ~ 1 .2 Mq, but clearly this is not the 
case with the '*'0/'^0 ratio. The purpose of this work is to try solving this problem, possibly also 
deriving further hints on the extra-mixing parameters that physical models must reproduce. 

In Sect. 2 we summarize the input data for our analysis, namely the spectra we used, the stellar 
parameters and the chemical analysis tools we adopted. Section 3 is then devoted to the description 
of the non-convective models assumed for explaining the newly determined C and O isotopic ratios 
of our two program stars. In Sect. 4 we then comment on the values found for the extra-mixing 
parameters and we derive on this basis some general conclusions. 



2. Input data 

2. 1 . Observed data and line lists 



Hinkle et al. 



(119951) . 



For Arcturus we used the electronic version of the Infrared Atlas Spectrum by 
We analyzed the ~ 5 yum-region spectra rationed to the telluric spectrum. Accurate wavelength po 
sitioning and identification of the main specie s contr ibuting in this region were recently performed 
in the solar infrared spectrum by 



Hase et al 



(I2010h . In particular, in the range 1800-2200 cm ' 
there are many weak and unblended '^C'^O and '^C'^O lines very sensitive to changes of the O 
isotopic ratios. In the case of Aldebaran, we used a spectrum in a similar but shorter spectral re- 
gion obtained on February 6, 1980 at the KPNO 4 m coude telescope using a Fourier transform 
spectrometer This spectrum was kindly provided by K. Hinkle. It has a spectral resolution of 0.016 
cm ', slightly lower than that of Arcturus (0.01 cm The spectrum of Aldebaran was cleaned by 



' The derivation of this ratio, mainly from optical bands (CN lines in the region ~ 8000 A) and from the 
infrared domain (CO lines at ~ 2.3 yum), is rather insensitive to the stellar parameters adopted. Indeed, these 
parameters affect almost equally the '^ '^CN and/or '^ '^CO lines. 

^ We report here only the O ratios derived by the above authors from the 5 /jm region, as in the 2 /jm region 

the available CO lines for K giants are usually weak and blended. 



20101) 



^ Among which the new measurement of '*0(p, a)'^N cross section provided bv lLa Cognata et al. 
These predictions do not change significantly with metallicity in the range -0.5<[Fe/H]< 0.0. (In the 

present work, we adopt the standard notation [X/H]= log(X/H)y,-log(X/H)0 where (X/H) is the abundance of 

the element X by number in the scale log (H)= 12.) 
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Table 2. Stellar parameters used to derive CNO abundances and isotopic compositions in a Boo 
and a Tau. 





a Boo 


a Tau 


Ref. 


Name 


Arcturus 


Aldebaran 




HR 


5340 


1457 




MK type 


K1.5 III 


K5 111 




reff[K] 


4290+ 50 


3981±75 


a 




196±21 


440±20 


b 


log (g [cms-^]) 


1.50±0.10 


1.20+0.30 


a 


R[Re] 


25.4±0.2 


45.2±0.7 


a 


[Fe/H] 


-0.50±0.07 


-0.13±0.13 


a 




1.7±0.1 


1.94±0.20 


a 


M[Me] 


1.08±0.06 


1.3±0.3 


b 



Notes. The f o llowin g ref erences are for Arcturus an d Aldebaran, respectively 



(a) 



Ramirez et al 



120121. 



1 2003) and 



Alves-Brito et al 



20101) 



(i) 



Ranurez & AUende Prietd d201lh : 



Rvdeetal 



l 2009h 



Lebzelter et al 



and references therein). 



telluric absorptions using the telluric spectrum of Arcturus (after some spectral resolution degrada- 
tion) with the IRAF task telluric. At the wave number position of the strongest telluric absorptions 
the removal was, however, unsatisfactory, thus these spectral regions were excluded in the analysis. 
A difficulty in this procedure is that many peaks are found above the level of unity in the rationed 
spectra. We detect these highest peaks (excluding those close to the regions with the strongest tel- 
lur ic abs o rption s) and fit a smooth curve passing through them following the similar method used 
by iTsuiil (120091) in Arcturus to place the continuum level in the rationed spectra. It is uncertain 
whether the continuum adopted in this way is a true continuum; however, different fits to these 
peaks resulted in very small differences (< 1%) in the continuum level. 

We made use of an improved molecular line lists in the ~ 5 jum region. Our list includes the 
molecules CO, C2, CN, OH, SiO, MgH SiS, H2O , being CO the main contributing molecule in 



the sp ectral region. CO lines come from lGoorvitchI (119941) : C2 lines are an update of lOuerci et al 



( 1197 lb (private communication); CN line s from B. Plez (pr ivate communic a tion), a s an update 



after the new energy levels c alculated by 



Ram et al 



are from 



Barber et al 



Skorv et al 



(I2OO6I) ; SiO ones from 



(120031) : OH ones from 



(l2010al) and 



Ram et al 



Langhoff & Bauschlichei 



Goldman et al 



(2010b): H2O lines 



1993 



Cami et al. 



; MgH on es from 



(120091) and OH 



( 1998 ). SiS onesfrom^ 
' 1 I 

ones from the HITRAN database dRothman & Gordonll2009D . The atomic lines are taken from the 



VALD v-0.4.4 database (iKupka et alEoOOh . A few line positions and intensities (mainly CO lines) 
were corrected b y comparing a theoretical spectrum of the Sun wit h the infrared solar spec trum 
( Hase etalJ20IO ). We used a MARCS atmosphere model for the Sun ( Gustafsson et allboOS ) with 
the solar abundances set from lAsplund et alj (120091) . For the C and O isotopic ratios in the Sun we 



also adopted the values suggested by these authors. The fit to the solar spectrum was excellent, in 
particular between 2100-2200 cm"' where most of the '^C'^O, '^C'^''^0 lines used to derive the C 
and O ratios were selected. 
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Table 3. CNO abundances and isotopic ratios 



Star 


los 


; 6(C) log 6(N) log 6 


-(O) '-c/"c '^o/'^o if-o/'^^o 


17Q/18Q 


a Boo 
a Tau 


8.06 ± 
8.25 ± 


0.09 (20) 7.67 ±0.13 8.76 4 
0.12(15) 8.05 ±0.11 8.48 ± 


0.17 9 ±2 (24) 3030±530(7) 1660±400 (18) 
0.14 10 ±2 (11) 1670±550 (6) 666±450 (9) 


0.55±0.12 
0.4±0.08 


Notes, log 


6(X) = 


log Hr/riH + 12, where log nx is 


the number density of element X. The number between 





parenthesis indicates the number of lines used. The N and O abundances are adopted from the literature. We 
exclude features with both a "O and "*0 contribution (see text). 



2.2. Atmosphere parameters 



a Tau and a Boo have been extensively studied with high resolution, high S/N spectra since 
1980. Their atmospheric parameters have been estimated by several authors using a variety of 
techniques. For Ajcturus, the s imple mean and stand ard deviation for the stellar parameters com- 



piled in the PASTEL database (iSoubiran et alJl2010l) is T^// = 4324 + 90 K, log g= 1.71 + 0.29, 
and [Fe/H]= -0.56 ±0.1, while for Aldebaran is Te/f = 3850 + 40 K, log g= 1.2 + 0.4, and 
[Fe/H]= -0.16 ± 0.1. Despite the considerable number of studies on these stars, the published 
parameters still do distribute randomly around the mean values, due to the impact of system- 
atic errors which vary among different studie^. We decided, therefore, to adopt the most recent 
determinations of the atmospheric parameters for both stars. They are based on high-quality vi- 
sual and infrared spectra and the use of MAR CS model atmosph eres (the same grid that we use 
here). For Arc t urus w e ad opted those derived in Rvde et al.l (120091) while for Aldebaran those from 



Ramirez et al 



( l2009l) . and 



Alves-Brito et al 



( 120 101) (see Table|2]l. The adopted values, in any case, 
do not differ significantly from the average values given in the R\STEL database. In the references 
quoted above the C, N, O, abundances were derived, as well as those for other species having an 
important role in the opacity of the model atmospheres (Si, Ca, Mg, S, Ti etc, see the original works 
for details). However, the atomic features present in the 5 fj.m region in both stars are very weak or 
severely blended; therefore we adopt the elemental abundances given by the quoted authors except 
for the CNO elements. We note that variations up to +0.25 dex in the metallicity of the stars has no 
impact in the derived C and O isotopic ratios. 

A spherical MARCS model atmosphere ( Gustafsson et alJl2008 ) was interpolated for each star 
from the grid of models for the parameters {Teff, log g, [Fe/H], ^mkro) in Table |2] We assumed 
for both stars 1 Mp and for Arcturus w e adopted an a-enhanced model {{afFt}- +0.4) as suited 



Avres & Linskv 



to its metallicity ([Peterson et alJll993l) . We note that Arcturus has a chromosphere (for example 



19751), so that the continuum flux in enhanced a t wavelengths shorter than about 



2000 A. This flux excess might be explained using a binary model dVerhoelst et al 



20051), but in this 



study we considered Arcturus as a single star, as the impact of the possible secondary companion 
is only important in the ultraviolet, a region that we did not use and that therefore does not affect 
our analysis. 
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2.3. Analysis of the C and O ratios 



For each of the model atmospheres adopted, synthetic LTE spectra were calculated in the region 
18 50-2200 cm^' with a s tep of 0.0002 cm ' using the TURBOSPECTRUM v9.02 code described 
in lAlvarez & Plea (119981) . and the hne list given above. The theoretical spectra were convolved 



with a Gaussian function with a FWHM~ 600 - 800 mA to mimic the spectral resolution plus the 
macroturbulence parameter. In order to estimate the C abundance and the '^C/'^C ratio, we selected 
a number of '^■'■'CO lines that are weak, unblended and, apparently, not affected by the procedure 
for removing telluric absorption. We note that variations in the N abundance by ~ +0.3 dex have 
no effect in the synthetic spectrum. Also, the '^C'^O lines are not very sensitive to changes to the 
O abundance by ~ +0 .2 de x, thus we decided to adopt in our stars the N and O abundances derived 

It is important to note that in the selection of 



by 



Rvdeetal 



( 120091) and 



Ramirez et al 



these CO lines we took into account the fundamental problem existin g in the 5 jjm spectrum of K 

Rvde et al.ll2002 ; Tsuil2009 ): namely, the fact that the CO 



Heaslev et aU1978 



giants (for example 

fundamental lines cannot be interpreted with a photospheric model only. These lines show an excess 
absorption (lines with equivalent widths log W/v > -4.75) which is probably non-photospheric in 
origin. Figures 1 and 2 show clearly that the central cor es of many CO line s cannot be reproduced 
using a ID photospheric model in a LTE calculation. iTsuiil (120081 120091) proposed instead that 
the extra-absorption originates from cool molecular layers referred to as a quasi-static molecular 
dissociation zone; this is sometimes named MOLsphere. Indeed, the formation of these molecular 
clouds in the outer atm osphere app ears to be a basic feature of all the red giant stars from early- 



K to late-M types (see iTsuii 



20091 for more details). Therefore, our estimates of the C absolute 



abundance and of the C and O isotopic ratios are based on a careful selection of the CO lines (i.e. 
log W/v < -4.75), so that in principle the approximation of LTE synthetic spectra using canonical 
ID photospheric models should be valid. By considering this, we derive in Arcturus a C abundance 
in agreement with th at in 



dex than the one by 



Rvde et alj (l2009b . In Aldebaran we derive a C abundance larger by 0.15 



Ramirez et al 



(l2009h (see Table O. On the other hand, the carbon isotopic 
ratios derived in both stars agree also with previous estimates, being 9 and 10 for Arcturus and 
Aldebaran, respectively (see Table O. 

Once the CNO abundances and the C isotopic ratio were estimated, in each star the absorption 
features due to '^O and ^^O were fitted by the synthetic spectrum varying the abundance of these 
isotopes in order to give the best fit to each feature, one at a time. We selected very carefully these 
lines to avoid as much as possible blending, cases where the position of the continuum was regarded 
as uncertain and/or the presence of weak telluric lines in the rationed s pectrum. This resul t ed in a 
lower number of useful '^O and '^O lines as compared to the study by [Harris & Lambert! (119841) . 



performed in the same spectral region. However, contrary to these authors, we did not assign any 
weight to any feature to compute the final O ratios. We excluded also features where both '^O and 
'^O were contributing. The abundances derived from the various features selected in this way were 
then combined to give a mean (Table 3). 



^ A recent di scussion on the curren t techniques to derive atmospheric parameters and the associated errors 
can be found in ] 



Lebzelter et al 



( 120121) . 



CN lines are very weak in the 5 yum region and cannot be used to derive the N abundance directly. 



7 



C. Abia et al.: Carbon and oxygen isotopic ratios in Arcturus and Aldebaran: 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T 




Q 2 I I I I I I I I I I I I I I I I I I I I I I I I I 

2140 2141 2142 2143 2144 2145 

Wavenumber (cm~') 

Fig. 1. Comparison of observed (dots) and synthesized spectra (lines) of Arcturus for different O 
ratios. Dotted line: '^c/'^c = 9 and '^O/^'O ="*0/"'0 = 0; continuos line: the same '^c/'^c 
with '^0/"'0 - 3030, '**0/"'0 - 1660. Note that the cores of the more intense CO Hnes are 
not reproduced by the theoretical spectrum (see text). Some of the C'^O and C'**0 lines used are 
marked. 

The major source of uncertainty in the derivation of the C and O ratios is the dispersion in the 
ratios obtained from the different Unes. Uncertainties due to errors in the atmospheric parameters 
(see Table |2]i are minor as compared to these. As mentioned before, changes in the metallicity of 
the model atmosphere up to +0.25 dex has no impact in the derived ratios. The same is true with 
changes in the N abundance by +0.3 dex and/or A Tp// - +100 K. Uncertainties in the gravity and 
microturbulence of the order of those quoted in Table|2] imply errors in the O ratios not exceeding 
about +100 in both stars. A larger impact on the final error comes from the uncertainty in the C 
and O abundances and the '^C/'^C ratio. All these sources of error, once added quadratic ally, give 
a total uncertainty of +180 for Arcturus, and +230 for Aldebaran. These figures may be safely 
applied to both O ratios. In Table[3]we indicate the total error in the C and O ratios after including 
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the dispersion in the ratios between different lines. It is evident that the dispersion from different 
lines accounts for most of the total error Systematic errors may be present, like for example the 
uncertainty in the continuum position and departures from LTE. Due to the weakness of the '^O 
and '**0 lines used, errors in the continuum position should affect the '^O and '**0 abundances 
almost equally, so that the ratio '^0/'**0 is probably more reliable. For the same reason, departures 
from LTE should be small in the layers where the key features are formed. 

When comparing our oxygen ratios with those derived in iHarris & Lambert! (119841) . we agree 
within the error bars only in the '^0/'**0 ratio in Aldebaran; we derive considerably larger 
'^O/'^O and '^O/^^O ratios in both stars. The differences in the atmosphere parameters adopted 
(Jeff, log g, ^) are not strong enough to explain the discrepancies. In fact by using atmo- 
sphere models from the Gustafss on et al. (2008) grid with the same stellar parameters and 
CNO abundances as adopted by Harris & Lambert (1984 ] we obtain: ^^0/"0^ 2325 and 
^^0/^0= 1430 for Arcturus, and ^^0/"0= 1540 and '^0/'**0= 560 for Aldebaran. The ra- 
tios are reduced significan tly but not enough ( note t hat we still agree in the '^O/^^Oratiojn 
Aldebaran). We recall that 



( 119761) and I.Tohnson et al 



Harris & Lamberd (11984) us ed atmosphere models from 



19801) . The 



Bell et al 



Bell et alj ( Il976l) 's models are in fact the ancestors of 



the new grid of spherical MARCS atmosphere models by Gustafsson et al. (2008). These new 
models considerably improve the atomic and (mainly) molecular opacity treatment, in partic- 
ular for giant stars, as well as many other physical approximations (see Gustafsson et al. 2008 
for a detailed discussion), so that we consider that new MARCS models mimic much better 
the real atmospheres of giants than the original ones by Bell et al. (1976) and Johnson et al. 
(1980). An additional source of discrepancy which is very probably the main cause of our dif- 
ferent findings is the line fist used for the CO molecule and the method employed for the con- 
tinuum placement. First, we note that lHarris & Lambert! (11984) only include CO fines in their 
synthetic spectra computations, while we consider a number of molecular and atomic species 
(see previous section) in the 4.5 /vm spectral region. Indeed, the CO molecule dominates 
the absorption in this spectral range, but we checked that the remaining molecular species 
(mainly CN, C2 and OH) introduce a veil of absorption which increases the line intensities and 
thus affects the isotopic ratios derived. In part icular the '^O/^^O and ^^O /^^O ratios increase. 
Secondly, the jgf-values of the CO lines used bv lHarris & Lambert! (!l984) were obtained from 



Chackerian & Tipping! (!l983h . while here we used those by Goorvitch (1994) on the basis of 



a more up-to-date electric dipole moment function. As a consequenc e, Goortvitch (1994) in 



Chackerian & Tip pins 



deed reported differences up to 3% in the A -values with respect to 
(!l983h and up to a factor ~ 6 in the dipole moments in some of the CO isotopes. On the other 
hand, note that the isotopic fines used are very weak so that the choice of the specific lines 
might affect the results, in a systematic manner (see above the discussion on the selection 
of the fines) if there is a systematic difference in the continuum placement. Definitely, the 



' These authors give only the C abundance derived ([C/H]= -0.7 and -0.3, respectively) without any 
indication of the N and O abundances derived/adopted nor error bars. Thus, in th is test we scale the N 
and O abundanc es according to the metallicity given in lHarris & Lambert! U984!) respect to the solar 



abundances from 



Lambert! 419781) (we infer that these solar abundances were adopted by these authors; 



note this circumstance in figures 5 and 7 below, where we adopted a conservative error bar of ±0.2 dex 
in the CNO abundances by Harris & Lambert) 



9 




Fig. 2. As Figure [U for Aldebaran (filled circles) in another spectral region. Dotted line: '^C/'^C 
= 9 and "O/'^O =i^0/"'0 = 0; continuous line: the same '^c/i^c with i''0/"'0 = 3030, i^0/"^0 
= 1660. Note again the difficulty to fit the cores of the strongest CO lines (see text). Some of the 
C'^O and C'**© lines used are marked. 



di scussion above might exp lain the differences between the O ratios derived here and those 



in 



Harris & Lambern ( 119841) . We beUeve that our figures are more reliable because they are 



based on the use of more accurate stellar parameters and CNO abundances, better atmo- 
sphere models and more complete and accurate line lists. 



2.4. Physical stellar parameters 

In order to check the reliability of the mass estimates available in the literature (see Table |2]l, we 
used the FRANEC code to construct our theoretical HR diagrams. (We warn the reader that the 
theoretical mass estimate for a single stellar object, obtained by fitting its HR diagram, is quite 
uncertain. On the other hand, more robust methods, such as the classical isochrones fitting, cannot 
be applied to our stars). This was done in the following way. i) First we selected the range of 
initial masses M,„, for which the theoretical HR diagrams could fit the L, Teff values measured 
and the available estimates for the stellar radii of our stars within their uncertainties, ii) Then, we 
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computed, over the mass range found, a grid of extra-mixing models, looking for the combination 
of Mini values and transport parameters that allowed the best global fit. 

With this procedure, we found very easy to reproduce the observations for a Tau, assuming 

initially solar abundance ratios (see section[3]); our best estimate for the mass agrees very well with 

previous determinations. For Arcturus, instead, things were considerably more complex. This is a 

slightly metal-poor star for which guessing the initial abundance ratios is not straightforward. For 

C and N we had to rely on the literature, as the measured data are certainly not the initial ones, 

being modified by FDU and extra-mixing. However, there are suggestions that Arcturus belongs 

to a peculiar stream of stars in the vicinity of the Sun; some of them date back to the early sev- 
J II 1 

enties (IEggenlll971h. Today this stream i s considered as being the relic either of an old dissolved 



cluster (|Ramfrez & 



dNavarro et al. 



Allende Prieto 



201 ll) or of a captured and disrupted dwarf spheroidal galaxy 



120041) . Its abundances seem to be rather similar to those of the Galactic thick disk, 
but cautions are mandatory. Moreover, there is a considerable dispersion of initial abundances in 
thick disc stars at the metalhcity of a Boo, so that we can only put weak constraints on them. 
The initial values adopted here include a C enhancement by -hO.2 dex and a N under-abundance 
by 0.1 dex ( Bensbv & Feltzin jl2006 : Matteucci & ChiappinJboOS ). With t hese choices, and using 
an initial enhancement of a-rich elements of -1-0.4 (including oxygen) from lPeterson et alJ d 19931) . 
the model HR diagrams that can fit the (L, Tejf) data within the uncertainties correspond to a mass 
interval from 1 to about 1.25 Mq. The final value we adopted (1.2 Mq) is the only one for which 
extra-mixing models can reproduce within the errors all the abundance information discussed in 
Section 1273] The uncertainties on the initial C, N data make however o ur solution less robust than 
for a Tau. Note that our value is slightly larger than the one derived by 



Ramirez & Allende Prieto 



2011 



(1.08 Mq), who used as reference the Yale code, with a lower a enhancement. 



3. Reproducing the abundances through an extra-mixing model 



Whate ver the mechanism is that drives non-convective mixing in red giants, IWasserburg et al 



(1 19951) showed that it can be approximated by a circulation occurring at a rate M, reaching down 
to a maximum temperature Tp, close to, but lower than, the H-burning shell temperature. In a 
diffusive approach the parameters would instead be the diffusion coef ficient D and the to tal mass 



involved: the two approaches can be shown to be roughly equivalent dNollett et al 



2003h in most 



cases. One has however to notice that, while this is certainly the case when the mixing speed is not 
relevant (so that the abundance changes depend only on a path integral of reaction rates), things 
might be different when the velocity of mixing becomes an important issue, as diffusive processes 
are always slow, while other transport mechanisms might not be so. In fact: i) the time available 
for mixing is not infinite; and ii) the nuclei involved are somet imes unstable with a re latively short 



half-life: se for example the case of ^Be, decaying into ^Li jPalmerini et al. 



2011b). On similar 



grounds, recently serious doubts have been advanced on some of the proposed mecha nisms, like 



rotation and thermohahne mixing dCharbonnel & Lagardel201C ; 



Palmerini et al. 



201 lb!) because of 



the small diffusion coefficients (or, alternatively, the slow mixing speeds) they can provide, which 
might make them inadequate to yield the observed abundance changes in the finite time assigned 



by the duration of the evolutionary stage (iDenissenkov & Merrvfieldll201 11) . 



1 
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Fig. 3. Left panel: a comparison of the observed L and T^jf values of a Tau with the adopted 
evolutionary track in the HR diagram, as derived from the FRANEC code. Right panel: zoom of 
the previous plot (top) and a comparison of model and observations for the L, R^mr relation. The 
main sequence (MS), first dredge-up (FDU) and (BLF) positions are marked. 



An important input to the models is the initial CNO admixture of the stellar composition. While 
for a Tau, a typical thin-disk red giant of relatively high metallicity, we can safely assume solar 
elemental ratios, as we already mentioned for a Boo this is certainly not the case (see discussion 
in Section l2!4l l. In our procedure, once the initial CNO abundances are selected, they are employed 
for choosing the opacity tables to be used. Enhancements in a elements and, to a lesser extent, in C 
easily introduce large changes in the models, from both the points of view of both nuclear physics 
(CNO burning efficiency) and radiative transfer (opacities). We underline this important point be- 
cause the theoretical HR diagram and the ensuing mass estimate strong ly depend on that. No te that 
using the most recent set of aZp/za-enhanced opacities (see for example Ferguson et al 12005 ) theo- 
retical curves are moved to redder regions of the HR diagram than previously found. Thus, we warn 
that all the mass estimates available so far in the literature are actually much more uncertain than 
curre ntly supposed (see for example 



Verhoelst et al. 



2005 



Tsuii 



2009; 



Ramirez & Allende Prieto 



201 ih . 

Our best solution implies, for a Tau M - I.SMq, [Fe/H] - -0.15; for a Boo, M - 1.2Mq, 
[Fe/H] - -0.5, [a/Fe]= 0.4, [C/Fe]= 0.2. As mentioned before, we performed a grid of calculations 
for both stars: for each choice of the mass allowing a fit (within the uncertainties) to the L, T^ff 
data, we derived the extra-mixing parameters trying to reproduce the chemical abundances; finally, 
we adopted the case allowing the best compromise in the fit of all the available data. There are no 
ambiguities on a Tau with this procedure, and the resulting HR diagram is shown in Figure[3] 

For a Boo, instead, the solution that we finally adopt appears to require a mass higher than so 
far assumed and is strongly dependent on the initial CNO assumed. This point is further commented 
later (Section [X2b . The results for the HR diagram and radius of a Boo are displayed in Figure 4. 

According to the above discussion, from the estimate of the mass we derive the time spent on 
the RGB from the bump in the luminosity function (the small dent in each RGB track indicated 
by the label "BLF") to the moment in which the observed values of L, Teff are attained. This is 
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Fig. 4. The same as Figure[3]for aBoo. Note the peculiar initial CNO content for this star (see text). 
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Fig. 5. The evolution of the '^C/'^C ratio and of the elemental CNO abundances as a function of 
the model effective temperature in selected extra-mixing runs performed for a Tau, as compared 
to observed data. Black st ars show the observa t ional d ata presented in this work, while open stars 
report abundances from [Harris & LambertI (Il984l) . A case without extra-mixing is reported for 
comparison (dotted line). See text for details. 



the time available for extra-mixing to operate: it turns out to be 46 Myr for a Tau and 37 Myr for 
a Boo. This knowledge allowed us to determine the mixing parameters on observational grounds 
(albeit with the mentioned cautions for a Boo). 



13 



C. Abia et al.: Carbon and oxygen isotopic ratios in Arcturus and Aldebaran: 




1500 



,1000 



500 



a Tau 


^ BLF 


: 1 





10 20 

12^/1 3c 



1500 



1000 



500 




500 1000 1500 2000 2500 



2500 
2000 
O1500 
°1000 
500 



30 









^ BLF: 


a Tau 









10 20 



30 



Fig. 6. Combinations of isotopic C and O ratios for selected extra-mixing runs performed for a Tau 
and for a case without extra-mixing (blue spot), as compared to observed data. The models and 
the meaning of the line types are the same as in Figure 5. The small squares indicate the envelope 
isotopic mix reached by each mixing case when the model effective temperature is comparable with 
the observed one. At the same mixing depth, the higher the mixing rate Me, the lower the reached 
value of the carbon isotopic ratio the low er is the value reached for t he carbon isotopic ratio. Open 
stars are the isotopic ratios derived by [Harris & LambertI ( 11984 . 



3.1. The technique of the computations 

In order to fix quantitatively the parameters of extra-mixing using as constraint s the o bservations 



NoUettetal 



(120031), and made 



of a Boo and a Tau discussed so far, we adopted the formalism by|] 
therefore parameterized calculations. After deducing the parameters that allow us to fit the mea- 
sured abundances, we derived the corresponding mixing speeds needed to achieve the observed 
abundances in the assigned time. On this basis, we could analyze which of the processes proposed 
so far in the literature offers a plausible physical mechanism for driving the mixing. 

In our procedure we adopted a post-processing code to compute the coupled phenomena of H- 
burning and transpo rt, taking the detailed stellar parameters from the output of the stellar evolution 
code (FRANEC: see 



; lCristallo et al .11201 In . which provides us with the physical structure of the star. 
For describing the nuclear physics phenomena coupled with dynamics one can simply use the 
total derivatives of stellar abundances: 



dNi__dNi_ dNj dM OR 
IF ~ 'dT'^ dM~dR~dt 



where the partial time derivative due to nucleosynthesis is: 
dNi 



dt 



+ NpNi-iA, 



(1) 



(2) 
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Fig. 7. The same as Figure 5 for a Boo. 



and the second term is due to mixing. Here the parameters A^, are the reaction rates and /l^/ 
are the decay ra t es. Fo r the nuclear par a meter s we adopted the recent upgrades suggested by 



Adelbergeretal 



(l201lh and 



Iliadis et al 



(l2010h . Details on our techn ique for compu t ing extra 



mixing and on the descending numerical computations are presented by 



Palmerini et al 



(12011 



3.2. The parametric results 



Leaving as free parameters the circulation rate M and the temperature Tp of the deepest layers 
reached by th e non-convective mixing , we could profit of previous work done for a wide sample of 
RGB stars bv IPalmerini et al.l (1201 1 alibi) for limiting the interval of variation of these parameters. 
Following the formalism used by the above authors, the circulation rate was expressed in units of 
10"^ Mo/yr, through the parameter M(,. Concerning the temperature of the deepest layers attained 
by the mass transport, Tp, this was considered through the parameter A - log Th - log Tp, Th being 
the temperature at which the maximum energy from H burning in the shell is released. Although 
this is certainly a not-very-intuiti ve way of expressin g the mixing depth, it offers a sort of rule- 
of-thumb criterion, established by 



NoUett et al 



(120031) . If A is lower than about 0.1, post-process 
mixing models are rather safe, in the sense that any nucleosynthesis occurring during the transport 
will not add significant energy to the stellar budget, thus not altering the reference stellar structure. 

On this basis Me was allowed to vary in the range from 0.01 to 0.3 and A in the range from 
0. 18 to 0.22. For each combination of values of the two parameters we computed the corresponding 
post-process mixing models, starting from the stellar structure along the RGB as provided by full 
calculations (with no extra-mixing) made with the FRANEC code. The technique adopted is that 
of reading from the outputs of the stellar code the physical parameters at and above the H burning 
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Fig. 8. The same as Figure 6 for a Boo. 
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Fig. 9. The same as Figure 4, but with a choice of 1 .08 Mq (see text for details). 



shell, up to the convective envelope, using them to compute the outcomes of the coupled processes 
of burning and mixing with the parameters adopted; the inputs from the full model stellar structure 
are refreshed periodically, in order to be sure to maintain coherence between the extra-mixing 
calculations and the real physics of the star. 

Among the many runs performed we show in Table 4: i) those computed for a Tau, using our 
best choice for its mass; and ii) those computed fo r a Boo, using our mass es t imate (1.2 Mq) and 



Ramfrez & Allende Prietol (1201 11) . 



iii) those run by adopting the mass values given by|] 

The comparison of model sequences for isotopic and elemental abundances with observed data 
is presented in Figures from 5 to 8 for the best cases selected from Table 4. Here the curves with 
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Fig. 10. Diffusion coefficients approximating our circulation models under the hypothesis of a dif- 
fusive treatment (see text for details). The parameters of the considered cases are shown in the 
labels. 

different Une types refer to different choices of the parameters, according to what is explained in 
the labels. 

In particular, Figure 5 shows the evolution of the '^C/'^C ratio and of the elemental CNO 
abundances along the RGB for extra-mixing models adopted for a Tau, as a function of the effective 
temperature (which in this case is a proxy for time, as red giants become progressively redder and 
cooler when they evolve). The label BLF identifies the abundances as determined by the FDU (these 
are the same as those characterizing the envelope layers at the Bump of the Luminosity Function). 
It is clear from the plot that there is no way of explaining the carbon isotopic ratio with a model 
that does not invoke extra-mixing. Such a case (dotted curve) presents in fact a too large ^^C/'^C 
ratio when compared to the observations. The same conclusion can be obtained by an inspection 
of Figure 6, where the model combinations of the C and O isotopic ratios are compared with the 
observations. In this plot the case without extra-mixing is represented by a single dot, since neither 
carbon nor oxygen change their isotopic ratios along the RGB. An extra-mixing model with a too 
deep penetration (those with A=0. 18) can be excluded on the basis of its large '^O/'^O ratios and of 
its low absolute carbon abundance (see Table|4|. A case with an intermediate depth (A=0.20) would 
alleviate the problem of oxygen isotopic ratios, but would still predict a too low carbon abundance. 
Thus, for a Tau a rather shallow extra-mixing is needed (A=0.22). Concerning its efficiency in 
terms of mass circulation, all the cases lie within the observed range, apart from that with the lowest 
efficiency (Mg = 0.01). The Mg - 0.012 case shows the maximum allowed '^C/'^C ratio, while 
extra-mixing models in the range 0.03 < M(, < 0.3 reproduce well all the observed constraints. We 
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remind again that the most critical point is certainly the '^C/'^C ratio, even if oxygen isotopic ratios 
help in constraining the extra-mixing depth . The choice of t he para meters found to be appropriate 
is in good agreement with the cases run by Palmerini et al. ( 201 lalfb ). so that a Tau appears to be 



really a typical template for Population I red giants, also from what concerns the mixing processes. 

Figures 7 and 8 integrate the previous picture by showing the more complex case of a Boo. As 
discussed previously, here the constraints from elemental abundances are rather weak, depending 
on assumptions made for their unknown initial CNO values. Despite this caution, it is noticeable 
that the extra-mixing parameters determined for a Tau (0.03 < Mg < 0.3, A = 0.22) also reproduce 
well the constraints for a Boo. This is, in our opinion, a very interesting finding, as the two red 
giants have a different metallicity. Many of the computed cases have been excluded on the basis 
of reasonings similar to those exposed for a Tau (see Table |4|i. Considering the full extension of 
the error bars for the observations, extra-mixing cases with a rather slow circulation (0.01 < M(, < 
0.03) cannot be a priori excluded for a Boo. Notwithstanding, the global quality of the fits is much 
better with larger Mg values, i.e. with a choice similar to that for a Tau. 

From a further inspection of Table |4] it emerges that, as far as the abundances are concerned, 
ad opting a slightly lower mass than found by us (for example for 1.08 Mq, the value determined 
by Ramirez & AUende PrietJl201 1 ) several (albeit not all) data can be reproduced with a different 
choice of the parameters (A ~ 0.18, M(, ~ 0.015), that is with a slightly deeper and slower extra- 
mixing. Remarkable discrepancies emerge in this case only for the absolute elemental abundances, 
whose modeling suffers for the uncertainties in the initial composition already discussed. However, 
in the case with M = 1.08 Mq a poorer fit to the luminosity and radius of a Boo is obtained (see 
Figure 9). We therefore maintain our previous choice of the mass (1 .2 M©) as our best case, but we 
remark that it, too, is rather uncertain. 



4. Discussion and conclusions 



We notice that the values of the extra-mixing parameters deduced for our stars nicely fit with 
those previously found to be typical for red giants of Population I, from their CNO isotopic and 

J — ~ — II — in 

elemental abundance ratios as well as from their Li content jPalmerini et al.ll201 1 alibi) . Under these 



premises, we are tempted to conclude that extra-mixing during the RGB phase presents common 
properties for solar or moderately low metallicities. The idea is suggestive, although we cannot 
at this level draw too firm conclusions. A wider sample of mod erately metal-poor, wel l measured 
red giants, permitting reliable statistics (like the one selected by Palmerini et al.ll201 la , for higher 
metallicities) would be required for this. 

It is relevant to compare our findings to those proposed i n the hteratu r e, ado pting 
the complementary view of a diffusive approach. In particular, jPenissenkovl (1201 Ol) and 



Denissenkov & Merrvfield 



(|201 ll) . recently performed ID studies of the effects of thermohaline 



diffusion (which can be easil y compared to our dis cussion) and then substantiated their results with 
2D- and 3D-simulations. In iDenissenkovl (120101) it was shown that thermohaline mixing might 
guarantee diffusion coefficients D„,u smaller than a few 10^ cm^/sec (see especially their Figure 
3, in which D„„i was normalized to the the thermal diffusivity. This last was assumed to be of the 
order of 10^, as specified in Table 1 of that paper). However, in the subsequent discussion (in par- 
ticular from Figures 7, 1 1 and 12) the author makes clear that, in order to fit all the red giant data 
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(including the values of C/Fe and the isotopic ratios '^C/'^C lower than about 15) one would need 
much larger values of Ahu - These last seem however to be compatible with thermohaline mixing 
only for low metallicities. In our case we can make a rough estimate of Dmix fro m the simpli 



Nollett et al 



fied co rrespondence between the circulation and diffusion treatments established by 
(bod), i.e. D,g,i, - (/ X M)/4npr\ 

For M = (0.3 - 3)x lO^^Mo/yr (the range found to be good for our stars) and adopting the values 
of the other parameters from the stellar code outputs we obtain values of Dequiv as those plotted in 
Figure 10 (see also the last column in Table |4] where average Df^„„, values of each extra-mixing 
case are reported). For the acceptable cases they cover a range centered around a few lO' cm^/sec 
(this last value being roughly the average). 

We can compare this with the data of Figure 12 in iDenissenkovl (120101) . where the curve for 



Dmix providing the best fit to the data is plotted as a function of the radius. The radius of the 
convective envelope border in a Tau, in the phases after the BLF, spans the range log[r/RQ] - 
-0.13 to -0.07. For these values the average value D„,u is again 10'' cm^/sec, in agreement with 
the value found by us. Hence, with a completely independent treatment, we confirm the results of 
the quoted paper. The diffusion coefficient required for explaining the observations of a Tau (and 
the more uncertain a Boo) must be quite large (Such a value for D„„ j implies in our cases a velocity 
of a few hundredths cm/sec). 

Note that also the treatment by|] 



Paknerini et al 



(1201 lah . favoring M values very similar to those 
of a Tau and a Boo, would provide the same consequences for the values of D„„ f . 

By pursuing a similar discuss ion in the framework of multi-D models of the radiative zones 



Denissenkov & Merrvfieldl (1201 ll) showed that this is related to the aspect ratios (length over di- 



ameter) of the dynamical instabilities generated in the simulations: with thermohaline mixing these 
aspect ratios would be too small, i.e. the unstable blobs would be too similar to "bubbles" instead 
of the required "finger"-like structures. 

Whatever approach is used, the results seem to converge in saying that pure thermohaline dif- 
fusion might have difficulties in explaining the observed abundances of high-metallicity red giants, 
at least when it is taken alone; t he po ssibihty of a modified magneto-thermohaline mixing was 

( 2009) but not yet substantiated by detailed models. When also 



Denissenkov et al. 



envisaged by 

the results bv lCharbonnel & Lagardd ( 12010 ') are considered, one sees that any diffusion induced by 
rotational effects is in its turn insufficient; in this case, actually, D is too small by several orders of 
magnitudes (see Figure 9 in that paper). 



We have to notice that in some recent works (see for example lAngelou et al. 



2011) the 



authors wisely avoid the use of the term "thermohaline diffusion" for indicating chemical 
readjustments started by a molecular weight inversion; they prefer the name Aju mixing. The 
use of the term diffusion could have remarkable impUcations. Indeed, while diffusion is intrin- 
sically a slow phenomenon and thermal diffusion should in fact be slow, it is not guaranteed 
that, in presence of an inversion of the /j gradient, the chemicals would diffuse with a speed 
comparable to that of heat. We have no elements at this stage to exclude that non-diffusive 
mixing of a suitable velocity might occur as a consequence of a "A//" effect: in that case, it 
would offer a realistic mechanism. Our analysis only underUnes that very slow mixing, as in 
diffusive processes, would be inadequate for explaining the chemical abundances; but there 
is clearly a lot we have still to learn about the real physical processes. 
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One ha s to u n derline that magnetic buoyancy, recently advoc ated in various ways by 



Busso et al. 



(120071) : 



Nordhaus et al 



(l2008h : 



Denissenkov et al 



(120091) is suitable to provide the 



mixing velocities (or the D,„u values) we req uire for the tw o stars examined. Looking back at old, 
seminal works by E.N. Parker (for example |Parkeiill974l) one sees that the velocity of buoyant 
magnetic structures, in presence of thermal exchanges with the environment, is roughly v = Kja^, 
where a is a typical linear dimension of the rising bubbles. Using parameters suitable for evolved 
red giants (Parker analyzed instead the case of the Sun), K turns out to be of the order of lO'"^ - 10'^ 
(cgs units). Hence, large structures (100-1000 km-size) would travel at moderate speeds, around 
the velocity required on the RG B, while small instab ilities (1 km-size) would provide the situation 
envisaged for the AGB phases in[ 



Busso et al 



(120071) . with high speed, close to the Alfven velocity. 
All the cases in between these extremes are possible. If this is the real physical situation (which 
fact has still to be proven on the basis of MHD simulations), then our results would suggest that on 
the RGB large magnetic domains, moving at moderate speed, are involved in the buoyancy. 

Note that, recently (IDrake et al.l 1201 ih . the association between non-convective mixing and 
magnetic activity seems to have been demonstrated n icely for the bright, very active RS CVn-type 
variable A And (see for example Andrews et aljl988 ). This star shows CNO anomalies well before 



reaching the BL R This finding, and the fact that one of our stars (a Boo) is known to have both 



a chromospher e (lAvres & Linskv 



1975|) and ind ications of photospheric magnetic fields from the 



Zeeman effect dSennhauser & Berdvuginall201 11) . seem actually to suggest magnetic mechanisms 
as very promising physical causes for driving extra-mixing in red giants. 
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Table 4. CNO abundances, C and O isotopic ratios and equivalent diffusive coefficients for the 
calculated extra-mixing models. 

Table 4 
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